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Using a data sample corresponding to an integrated luminosity of 342 fb -1 collected with the 
BABAR detector at the SLAC PEP-II electron-positron storage ring operating at a center-of-mass 
energy near 10.58 GeV, we measure B(r" -> -k~ -k+ v T ) = (8.83 ± 0.01 ± 0.13)%, B{r~ -> 
K-TT'n+Ur) = (0.273 ± 0.002 ± 0.009)%, B{t~ K~-k-K+v t ) = (0.1346 ± 0.0010 ± 0.0036)% 
and B(t~ -> K~K~K + v T ) = (1.58 ± 0.13 ± 0.12) x 10~ 5 , where the uncertainties are statisti- 
cal and systematic, respectively. Events where the n + n~ pair is consistent with coming from a 
Kg are excluded. These include significant improvements over previous measurements and a first 
measurement of B(t~ — > K~ K~ K + v T ) in which no resonance structure is assumed. We also re- 
port a first measurement of B{t~ —* 4>n~v T ) = (3.42 ± 0.55 ± 0.25) x 10 , a new measurement of 
B(t~ -^(f>K~v T ) = (3.39±0.20±0.28)xl0 -5 and a first upper limit on B(t~ -» K~ K~ K + v r \ex..tj>]). 

PACS numbers: 13.35.Dx, 14.60.Fg, 14.40.Cs, 14.40.Ev, 12.15.Hh. 



The weak interaction coupling strength between the 
first and second quark generations [l[ can be probed with 
unprecendented precision in hadronic t decays having 
net strangeness of unity in the final state using e + e~ — > 
t + t~ data collected at the e + e~ B-factories. Inclusive 
measurements of the strange spectral function, obtained 
from t lepton decays to final states containing kaons, 
provide a direct determination of the strange quark mass 
and Cabbibo-Kobayashi-Maskawa (CKM) mixing ele- 
ment \V U s\ 0- A significant improvement on the pre- 
cision of B(t~ — > K~~k~ti^v t ) in particular, where one pie of t~ 
measurement Q disagrees by more than two standard de- 
viations from the others [US]) will have the most immedi- 
ate impact on the determination of these two fundamen- 
tal Standard Model (SM) parameters using r decays @. 
Measurements of t~ — > <f>Tt~v T and — > <pK~y T pro- 
vide an interesting laboratory for studying Okubo-Zweig- 
Iizuka (OZI) suppression (in- 
significant improvements on measurements of B(t~ — > 
7T — 7r _ 7P ^ T ), B(t~ — > K~ti~it + i' t ) and B(t~ — > 
K~tt~ K + is T ) are reported together with a first measure- 
ment of B(t~ — > K~ K~ K + v T ) without resonance as- 
sumptions (charge conjugate decays are implied). The 
data sample corresponds to an integrated luminosity of 
L = 342 fb -1 recorded at a center-of-mass (cm.) en- 
ergy (y/s) near 10.58 GeV using the BABAR detector 
at the SLAC PEP-II asymmetric-energy e + e~ storage 
ring. With a luminosity-weighted average cross-section 
of cr e+e -_ +T+T - = (0.919 ± 0.003) nb @, 0, this corre- 
sponds to the production of 3.16 x 10 8 r-pair events. 

The BaBAR detector [Io[ has a silicon vertex tracker 
(SVT), drift chamber (DCH), ring- imaging Cherenkov 
detector (DIRC) and electromagnetic calorimeter (EMC) 
all contained in a 1.5-T solenoid. The iron flux return of 
the solenoid is instrumented to identify muons. 

r-pair events are simulated with higher-order radiative 
corrections using the KK2f Monte Carlo (MC ) g enera- 
tor with t decays simulated with Tauola [i"llfl2| using 
measured rates [13j. The detector response is simulated 



with GEANT4 14J. Simulated events for signal as well as 
SM background processes d, 11, [H, 15, are recon- 
structed in the same manner as data. The MC samples 
are used for selection optimization and systematic error 
studies. The number of simulated non-signal events is 
comparable to the number expected in the data, with 
the exception of Bhabha and two-photon events, which 
are not simulated. Data studies demonstrate that these 
backgrounds are negligible. 

The basic analysis strategy is to select a pure sam- 
• h~h~ h + v T decays from e + e~ — > t + t~ 
events by requiring the partner r + to decay leptoni- 
cally. Within this sample, each of the /i mesons is 
uniquely identified as a charged pion or kaon and the de- 
cay categorized as t~ — > 7r~7r~7r + i; r , t~ — > K~ tt~ ir + v T , 
t~ — > K~ir~ K + v T or t~ — » K~ K~ K + v T . An effi- 
ciency migration matrix, £ t j , initially obtained from MC 
simulations, is used to correct for efficiency losses from 
all stages of event selection and includes cross-feed be- 
tween the four signal channels where i (J) is the se- 
lected (true) decay mode index. The £y matrix is modi- 
fied using data control samples of kaons and pions from 
D* + — ► tt + D°, D° — ► tt + K~ decays to account for small 
differences between MC and data. The number of decay 
mode j signal events measured in the sample, N^ lff , is 
then: 



N 



Sig _ 



-1\ / lyr Data 



N 



Bkg 



(1) 



where Nf ata is the number of data events selected in 
decay channel i and N^*" 9 is the estimated number of 
background events in decay channel i arising from sources 
other than t~ — > h~h~hrv T . The branching fraction for 

decay mode j is then Bj — 2Ca - 3 — . 

The r~ — * lrrlrrh + v T sample is obtained by selecting 
events with four well reconstructed tracks having zero to- 
tal charge, where none of the tracks originate from the 
conversion of photons in the material of the detector. All 
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four tracks are required to lie within the geometrical ac- 
ceptance of the EMC and DIRC and, to ensure the tracks 
can reach the DIRC, are required to have a transverse 
momentum of at least 250 MeV/c. If any two photons in 
the event are identified as coming from 7r° — * 77, the 
event is removed from the sample. The event is divided 
into hemispheres in the cm. by a plane perpendicular to 
the thrust axis [13] ■ One of the hemispheres, the 'lepton- 
side', is required to contain a r — > evv or r — > [ivv decay. 
If the hemisphere with an electron (muon) has any neu- 
tral particles with more than 1.0 (0.5) GeV the event is 
not selected. The e + e" — > qq backgrounds are further 
reduced to ~ 0.1% of the r~ — ► h~h~h + v T sample by 
requiring a thrust magnitude > 0.85. 



Backgrounds from Bhabha and e H 



l-i ji events 



with photon conversions are suppressed by requiring the 
momentum of the lepton-side track to be less than 80% 
of \/s/2 and a requirement that none of the three signal- 
side tracks pass an electron identification algorithm. Re- 
maining non-r background events, including those from 
two-photon processes, are reduced by requiring the event 
missing cm. transverse momentum to be > 0.009-\/s. 

The remaining background in the sample is predom- 
inantly from other t decays containing ir° and Kg 
mesons. Events containing a K s are identified and re- 
moved. Residual backgrounds from decays having a tt° 
are reduced by requiring the total energy in the signal- 
side hemisphere deposited in the EMC which is unas- 
sociated with the three charged hadron tracks to be 
<200MeV. With these requirements, contributions from 
r~ — > h~niT v and r~ — > hrhrh^lifiv are negligible. 

A track in the r~ — > h~h~h + v T sample is classified 
as a kaon using a likelihood approach to combine infor- 
mation from the DIRC, DCH and SVT with a charac- 
teristic kaon identification efficiency of ~ 80% and pion 
misidentification rate of ~ 1%. If the track fails to be 
identified as a kaon, it is classified as a pion. Events are 
selected if the signal-side decays are identified as t~ — > 
7T _ 7T — 7P~i/ T , t~ — * K ~7r~7r + i/ T , t~ — > K~ ir~ K + v T or 
t~ — > K~ K~ K + v T with decays having a wrong charge 
combination removed. 

The diagonal elements of £y excluding particle iden- 
tification (e), numbers of selected, background and sig- 
nal events determined using Eq. [T] are shown in Ta- 
ble U for all modes. The increase in e with number 
of kaons in the decay is associated primarily with the 
transverse momentum requirement and the Kg veto. 
The background fraction from r decays to 7r~7r~7r + 7r°i/ 
(K-ir-ir+ir a v, K'Tr' K+n°v, K'K-K+^u) in the 
7r~7r~7r + (K~7r~ir + , K~ir~ K + , K~K~K + ) candidate 
sample is estimated to be (3.6 ± 0.3)% ((2.3 ± 0.4)%, 
(0.4 ± 0.1)%, < 5.0%). Non-r backgrounds comprise less 
than 0.5% of each channel's final event sample. 

The component of £y associated with the particle iden- 
tification, My, is shown in the first four rows of Table HT1 
Note that this matrix includes efficiency losses associated 



with cross-feed of the wrong charge combinations and 
small factors associated with data control sample cor- 
rections to the MC cross-feed efficiencies, therefore the 
columns of the table are not expected to sum to 100%. 

Systematic uncertainties are assigned for: luminosity; 
cross-section; migration matrix elements, which includes 
MC statistical and systematic errors associated with the 
efficiency and particle identification; signal-mode model- 
ing; modeling of the EMC and tracking response includ- 
ing scale and resolution uncertainties, the sensitivities 
of the measurements to the modeling of hadronic and 
electromagnetic showers in the EMC, and tracking effi- 
ciency; modeling of the trigger; and modeling of the back- 
grounds, including uncertainties on cross-sections and 
branching fractions. These are summarized in Table IIIII 
along with the total error correlation matrix. The abso- 
lute normalization is a significant source of the correla- 
tions. 

The branching fraction results of this analysis are 
presented in Table U together with the world aver- 
age values/limit published by the Particle Data Group 
(PDG) [13j, with which they are consistent and sig- 
nificantly more precise. In all four channels, the re- 
sults where the lepton-side has an electron are consistent 
with those where it has a muon. Our measurement of 
B(t~ — > 7T _ 7r~7r + z/ r ) is also consistent with a precision 
B(t~ — > 7r _ 7r _ 7r + [ex.a;]) measurement 18( after account- 
ing for the u>. Our measurement of B(t~ — ► K~tt^tt + u t ) 
is in agreement with Q and disagrees by more than two 
standard deviations from 0] and [Bj]. We report a first 
measurement of B(t~ —* K~ K~ K + v T ) in which no res- 
onance structure is assumed and which has a significance 
in excess of 8er. 

A (/>(1020) contribution is seen in both the t~ — > 
K~n~ K + v T and r~ — > K~ K~ K + v T decay modes. The 
use of a kaon selection algorithm with higher efficiency, 
but less purity, provides significantly higher signal-to- 
background for r~ — > <f)K~v T . The t~ — > <f>Tt~v T 
{t~ — > <f)K~v T ) signal has a 5.7a (9.8a) level of signifi- 
cance. The last two rows of Table [TT1 lists the My matrix 
for this higher efficiency selection. The K + K~ invari- 
ant mass distributions for the t~ — > K~n~ K + v T (t~ — > 
K~ K~ K + v T ) mode using this sample is shown in Fig.Q] 
©. Below 1.09 GeV (1.15 GeV), after background sub- 
traction of the iioyi-K + K~ events, the K + K~ invariant 
mass distribution from the r~ — > 4>Tt~v T (r~ — * <j}K^u T ) 
decay is well described by a Breit-Wigner function convo- 
luted with a Gaussian resolution function for the signal 
and a third-order polynomial (function in fl9| ) for the 
background and is used to fit for the number of events. 
A binned maximum likelihood fit yields 344±42 (274±16) 
t~ — > 4>ir~v T (r~ — > <f)K~v T ) candidates. MC estimates 
of the subtracted qq background events with a <j) con- 
tribute a 5.2% (0.7%) uncertainty. The 4>tt~ candidate 
sample has an additional 4.3% uncertainty arising from 



potentially peaking r 



background. The 
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TABLE I: Diagonal elements of £ij excluding particle identification, numbers of selected, background and signal events deter- 
mined using Eq. [T] The branching fraction results for this work are presented in the fifth line where the first error is statistical 
and the second systematic. The PDG average branching fractions/limit fl3j ] are presented in the sixth line. 





T — > TV TT TT + U T 


T — > K IX TT + V T 






e 

jyBkg 
J\[Sig 

B (this work) 
B (PDG average) 


0.028 

(f .5953 ± 0.0013) x 10 6 
(0.0642 ± 0.0002) x 10 6 
(55.59 ±0.05) x 10 6 
(8.83 ±0.01 ±0.13)% 
(9.13 ±0.05 ±0.46)% [4] 


0.031 

(6.956 ± 0.0026) x 10 4 
(0.2263 ± 0.0064) x 10 4 

(171.5 ± 1.2) x 10 4 
(0.273 ±0.002 ±0.009)% 
(0.33 ± 0.05)% 


0.035 

(1.819 ±0.013) x 10 4 
(0.0145 ± 0.0008) x 10 4 
(84.71 ± 0.66) x 10 4 
(0.1346 ± 0.0010 ± 0.0036)% 
(0.154 ±0.009)% 


0.039 
275 ± 17 
2.5 ± 1.5 
(9.93 ± 0.84) x 10 3 
(1.58 ±0.13 ±0.12) x 10~ 5 
< 3.7 x 10- 5 @90%CL 



TABLE 11: Component of the efficiency migration matrix as- 
sociated with pion-kaon particle identification, Mij, in per- 
cent. The t~ — > 4>tt~v t and r~ — > <f)R~v T study employs a 
higher efficiency /lower purity kaon selection algorithm. 



Candidates 


Decay modes 
TTTvn Rmr RttR RRR 


TVKTT 

Rnir 

RttR 

RRR 


97.40 22.49 4.73 1.02 
1.42 74.87 16.43 6.38 
0.01 0.49 59.63 25.54 
- 0.26 50.87 


(f>TT 

cj,R 


72.54 19.20 
0.83 66.06 



TABLE III: Upper: Systematic uncertainties (%). Lower: 
Correlation matrix from stat. syst. covariance matrix. 





7V7V7T 


Rtvtv 


RttR 


RRR 




0.9 


0.9 


0.9 


0.9 


e~ — >r + 1 — 


0.3 


0.3 


0.3 


0.3 


Mij and particle ID 


0.4 


3.0 


1.9 


4.9 


signal modeling 


0.2 


0.2 


1.3 


2.0 


EMC and DCH response 


0.8 


0.9 


0.8 


1.2 


trigger 


0.1 


0.1 


0.1 


0.1 


backgrounds 


0.4 


0.7 


0.4 


5.5 


Total 


1.4 


3.4 


2.7 


7.8 


TVTV7V 




0.544 


0.390 


0.031 


Rtvtt 






0.177 


0.093 


RttR 








0.087 



fit parameterization contributes a 1.0% (2.0%) error. Ac- 
counting for B{(j> K+K~) [H], B(r~ -► <j)iT-v T ) = 
(3.42 ± 0.55 ± 0.25) x 10~ 5 and B{t~ -► §K~v T ) = 
(3.39 ± 0.20 ± 0.28) x 10~ 5 with a correlation of -0.07. 
From the fit we find no evidence for t~ — * K~ K~ K + v T 
without a <j> and set a first upper limit on B(t~ — » 
K-K-K+v T [ex.(p]) < 2.5 x 10" 6 at 90% CL. 

This is the first measurement of B(t~ —> <^tt^v t ). It 
is consistent with a CLEO limit and larger than the 
(1.20 ± 0.28) x 10 -5 value predicted by a vector meson 
dominance model Q- Our B{t~ — > <pK~v T ) measure- 




FIG. 1: R + R~ invariant mass in -> R~ tt~ R + v T decays. 
Data are represented by points with error bars, MC of r~ — * 
R~tt~ R + v T by the open histogram, cross-feed from other 
r~ — > h~h~h + v T channels by the light shaded histogram, 
and non-r backgrounds by the dark shaded histogram. The 
inset shows the background-subtracted data with the fit (solid 
line) and non-resonant component (dashed line). 

ment is consistent with a recent Belle result [2ll. Recent 
calculations using a meson dominance model [22j agree 
with our B(t~ — > <pir~v T ) and B(t~ — > <pK~v T ) mea- 
surements and with the ratio §^f=^fic=£\ — 0.99±0.21. 
Our measurements of B{t~ — » K ~ tt~ tt + v T ) and B(t~ — > 
K~ K~ K + v T ), when combined with improved measure- 
ments of the other strange decays, will constrain the 
CKM element | V us | better than unitarity bounds @, [23| . 
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FIG. 2: K + K~ invariant mass in r~ — > K~ K~ K + v T de- 
cays with two entries per event. Data are represented by 
points with error bars, the open histogram is the MC of 
r~ — > cj>K~v T , the light shaded histogram is the cross- 
feed from the other r~ — > h~h~h + v T channels, primarily 
from r~ — > K~ tt~ K + v T , and non-r backgrounds by the 
dark shaded histogram. The inset shows the background- 
subtracted data with the fit (solid line) and non-resonant 
component (dashed line). 
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